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It is difficult to describe in a few pages the numerous specific techniques 
used to study absorption lines seen in QSO spectra and to review even rapidly 
the field of research based on their observation and analysis. What follows is 
therefore a pale introduction to the invaluable contribution of these studies to 
our knowledge of the gaseous component of the Universe and its cosmological 
evolution. A rich bibliography is given which, although not complete, will be 
hopefully useful for further investigations. Emphasis will be laid on the impact 
of this field on the question of the formation and evolution of galaxies. 

1. INTRODUCTION 

Although steady progress is made towards detecting faint objects in emission, 
the high-redshift objects detected this way up to now are drawn from a partic- 
ular population of powerful emitters (e.g. |8S|). On the contrary, absorption 
may reveal standard objects such as a normal galaxy or intergalactic gaseous 
clouds close enough to the line of sight to QSOs all over the range of accessible 
redshifts from zero to the highest QSO redshift z — 4.897 j78|. A typical spec- 
trum is shown in Fig. [j]. Absorption line systems observed in QSO spectra are 
generally classified into three categories: 

(1) The metal-line systems in which a large number of elements, in different 
ionization stages, is observed, from C I to C iv or O I to O VI. The Lyman 
limit systems (LLS) are optically thick at the H I Lyman limit (912 A) and 
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Fig. 1. — QSO spectrum from Songaila |5o| obtained with the high-resolution spec- 
trograph on the Keck telescope. The emission redshift of the quasar Q1159+123 is 
z cm = 3.502; its strongest emission lines are Ly/31025/Ovil036 at 4615 A, Lyal215 
at 5473 A, Nvl240 at 5582 A, Siivl400 at 6300 A and Civl549 at 6974 A. To the 
blue of the quasar Lya emission line, the number of absorption lines is very large. 
Most of the lines are attributed to Lyal215 from neutral hydrogen in intergalactic 
clouds; when the quasar is at large redshift like this one, the corresponding Ly/3 lines 
are redshifted in the optical band and observed, confirming the identification. The 
break in the continuum at 4150 A reveals the presence of a cloud optically thick 
at the Lyman limit (912 A at rest). To the red of the quasar Lya emission line, 
the number of absorption lines is smaller. These lines can be attributed to different 
ionization stages of metals. Since most of the strongest absorption lines are in the 
near UV at rest, a wide range of ionization stages is accessible for the analysis of the 
physical conditions in the gas. 



have therefore H I column densities iV(H i) > 2xl0 17 cm~ 2 . They are closely 
related, at any redshift, to galaxies. The latter have been detected by direct 
imaging and follow-up spectroscopy at low and intermediate redshift |8(| . 

Studying their number density and physical properties (kinematics, ionization 
state, abundances) is then a unique tool to trace galaxy formation. Among 
these systems, the rare damped systems, characterized by a very large H I 
column density (A^(H i) > 2xl0 20 cm -2 ) are often considered to be associated 
with galactic disks or protogalactic disks p8|. 

(2) The Lya lines with no metal lines detected at the same redshift, are of 
intergalactic origin at high redshift (e.g. ]76fl) but could somehow be associated 
with galaxies at low redshift (4^]. Information about these systems (column 
density, Doppler parameter, kinematics, clustering) is derived from line-profile 
fitting. The latter requires high-resolution high-quality spectra that are difficult 
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to obtain on 4m-class telescopes. Recent studies of absorptions coincident in 
redshift in the spectra of gravitationally lensed QSO images Q or QSO pairs 
(e.g. |22|]) have shown that the mean radius of the absorbers is surprizingly 
large (up to 500 kpc). The idea that the Lya gas could trace the potential 
wells of dark matter |l0| and reveal its characteristical filamentary structure 
has been successfully investigated using iV-body simulations |5g, |64j . 
(3) The broad absorption line (BAL) systems are characterized by impressive 
absorption troughs from different ions of low and high excitation, extending 
from up to 60000 km s _1 outflow velocity relative to the QSO emission 
redshift. It is widely accepted that the gas is very close to the centre of the 
QSO host-galaxy and may be part of the broad emission line region (|9(| |92]| , 
see also [^3| ). They thus are intimately related to the AGN phenomenon. 
Although AGN activity is an important factor of galaxy formation, the role of 
BAL systems in galaxy formation is unclear at the moment. This is the bad 
reason why they are not described in this lecture. 



2. ABSORPTION LINES 



In this Section, we summarize a few basic and classical results on the formation 
and characteristics of absorption lines, that are important to keep in mind when 
discussing the scientific issues presented in the following Sections. 



2.1. Definitions 

An absorption line has an equivalent width w bs in A defined as: 

^—-!-d\= I "(1 - e- T ^)d\ (1) 



W Q bs 

where I is the observed spectral intensity, I c the interpolation of the absorption- 
free continuum over the absorption feature and r(A) the optical depth. It is 
apparent from Eq.(Q) that the observed equivalent width does not depend on 
the spectral resolution. For redshifted absorption lines, u> b s = Wvest x (1+^abs)- 
If we assume that the absorbing atoms have a gaussian velocity distribution of 
mean v a relative to the referential R, then it is possible to estimate the optical 
depth in Eq.(Q) as seen by an observer in R. An atom with velocity v absorbs 
a photon of frequency v in R with a cross-section a(v') where v' = v/(X-v/c), 
v taken to be positive for an atom that goes away from the observer and, 

1 f + °° , „ <— >> 2 



N-=r / cr(z/)e &~ dv (2) 
v 7 ™ 



DC 



where N is the column density, the total number of atoms per surface unity 
integrated over the gaussian distribution centered at v Q on the line of sight; b 
is the Doppler parameter, b = v Tms V2 = FWHM/2VLn2; FWHM is the full 
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width at half maximum. When the velocity field is not turbulent, the Doppler 
parameter is related to the temperature of the gas by 



V=^ = 12.8(kms-^§ (3 ) 

where m is the mass of an atom, M its mass number and T 4 = T/10 4 K. In 
case the turbulence can be modelled as a gaussian velocity distribution with 

quadratic mean velocity y§V and mean velocity 2V/ v^, then b 2 ot — b 2 h + V 2 . 
The cross-section for the absorption in an atomic transition is characterized 
by an absolute (classical) value, modulated by the oscillator strength / and a 
frequency dependence which is a consequence of the finite life-time of the upper 
level: 2 7 

° = f x ~A x -1 — T-TT^ ( 4 ) 

47re m c c tt {v - v ) 2 + (^) 2 

where the usual notations for the fundamental physical constants have been 
used and 7 is the total damping constant (de-excitation rate of the upper 
level). Eq.(||) reduces to 

r(A) = 1.498 x l(T 2 ^^ff(a,u) (5) 

with 

rr/ \ a f + °° e~y 2 , A 7 c f , v. A \ 

H(a,u) = - ^—^dy; a=-\; u = -- 1 + - -- 6 

nJ-00 {u~yY + a 2 Anb b\ c A Q J 

H(a, u) is called the Voigt function which is a convolution of a gaussian function 

and a lorentzian function. 

The optical depth at the center of the line is: 

r.- 1.497 xlB- ^™-'^ (7, 
o(km s" 1 ) 

More details can be found in e.g. [B6j. 



2.2. The Voigt Profile 

Under the assumption that the velocity distribution of the atoms is described by 
a Gaussian function, the overall shape of the optical depth r over the absorption 
line is a Voigt function. The core of the Voigt function is gaussian and the far 
wings are lorentzian. The absorption feature itself, which is simply e _T , is 
called a Voigt profile. The different regimes are illustrated on Fig. 
The absorption of photons with wavelength Al, or on a velocity scale Vi = c(Ai — 
Ao)/Ao, is due to a continuous sum of contributions by the atoms with different 
velocities. The probability that atoms 1 of velocity absorb the photons Ai 
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is given by the value of the lorentzian at V±. Since Vi is close to the center of 
the absorption line, it is apparent that atoms 2 of velocity V± will dominate 
the absorption because the variations of the lorentzian is much faster than the 
variations of the gaussian from to Vi . The variations of the optical depth 
in the core of the line will thus be governed by the gaussian function. In con- 
trast, the absorption of photons with wavelength A2, or on a velocity scale 
V<i = c(A 2 — Ao)/Ao, is dominated by atoms with velocity 0. Indeed, far from 
the center of the line, the variations of the gaussian function are much faster 
than the variations of the lorentzian. The number of atoms with velocity V2 is 
negligible. If absorption is seen at A2, this is due to the non-negligible prob- 
ability of an atom 1 to absorb a photon with wavelength A2. Of course, the 
total column density must be large for the optical depth to be non negligible 
in the wing. Thus when using the sketch of Fig. ||, it must be noticed that 
the transition between both regimes is very sharp around u ~ ±3.2 (see Eq.^) 
which corresponds to more than twice the FWHM. This is why the damping 
wings are prominent only for strong lines. 

Of course the assumption of a gaussian velocity distribution can be questioned. 
Under other assumptions, the absorption profile can be significantly different 
from a Voigt profile. This is particularly important when discussing the char- 
acteristics of the absorptions derived from line profile fitting (see |4l], 47 p8|). 



2.3. The Curve of Growth 

From Eqs.(||)(^)(|l]) one can calculate the curve of growth which is the relation 
between the equivalent width of a line and the column density for different 
values of the Doppler parameter. The curve of growth for the H 1 Lya transition 
is plotted on Fig. | for b = 5, 10, 20 and 30 km s _1 . 
There are three distinct regimes: 

• When the column density is small (r Q < 0.1), the absorption line is optically 
thin and the equivalent width does not depend on b. This is the linear part of 
the curve of growth, where the determination of TV from w is easy and reliable. 
For any transition, 

N(cm- 2 ) = 1.13 x 10 20 ^^- (8) 

Fig. H shows that for H IA1215, a log N = 13 line is partially saturated whatever 
the b value is. 

• The logarithmic or flat-part of the curve of growth is characterized by the 
large dependence of N on b at a given w. In this regime, the determination 
of b and N are very uncertain except when several lines of the same ion are 
used. Equivalent width and optical depth at the center of the line (see Eq.|?]) 
are related by: 



— = 2-y/Ln(r ) (9) 
A c 

• The absorption lines that are on the saturated part of the curve of growth 
are characterized by prominent damping wings. The equivalent width does not 
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Absorption at VI 

1: gauss(0)*lorentz(Vl) 
2: gauss(Vl)*lorentz(0) 



Absorption at V2 

1: gauss(0)*lorentz(V2) 
2: gauss(V2)*lorentz(0) 



Atoms 2 dominate -> gaussian shape 



Atoms 1 dominate -> lorentzian shape 



Fig. 2. — Illustration of the formation of the Voigt profile: over the profile, the optical 
depth has a gaussian shape in the core of the absorption line and lorentzian shape 
in the extended wings (see the text). The functions have arbitrary widths and are 
normalized to one at the maximum. 
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Fig. 3. — Curve of growth: logarithm of the equivalent width (to in A) versus loga- 
rithm of the column density (N in cm~ 2 ) for different values of the Doppler parameter 
o in km s _1 ). The curves are calculated for H I Lyal215. The three characteristic 
regimes are illustrated (see the text). 
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depend on b and the column density determination from measurement of w or 
fitting the wings is very accurate. In that case (see Eqs.^ and [7]), 

f = 2.64^ x V~r (10) 

For H IA1215, 7 = 6.265xl0 8 s" 1 and: 

iV(cm~ 2 ) = 1.88 x 10 18 w r 2 cst (A) (11) 

Fig. H shows clearly that, for realistic values of b, the H IA1215 transition is 
in this regime for log N > 19. The definition of damped systems as absorbers 
with log iV(H 1) > 20.3 is artificial and was introduced to search for damped 
candidates in low spectral resolution data (see |J7|). Indeed, for such column 
densities, the absorption line has w b s > 15 A at z > 2. The probability that 
such a feature be the result of blending of weaker lines is small. For smaller 
column densities, the probability for such confusion is much higher. 
The values of A and / are given for illustration for a few transitions in Table 1. 
A more complete list can be found in |57]] |Q |77j . 

Note that most of the transitions are from the true ground state. Exceptions are 
found (i) when the excited levels of the ground-state are populated by pumping 
from the true ground-state by the background radiation (C i*AA1560,1656, e.g. 



82 , this leads to determinations of the CMB temperature at high z; C n*A1334, 
e.g. [Q, is used to derive limits on the electronic density) and (ii) in a few 
associated or BAL absorption systems in which the density is high enough for 
upper levels to be populated by collisional excitation (e.g. pq| ). 

The column density and Doppler parameter for an ion can be formaly deter- 
mined if two transitions from the same lower level are observed on the flat part 
of the curve of growth. Indeed from Eqs. (|9j) and (0) it is easy to extract N and 
b. This is the doublet method which has been extensively used because dou- 
blets are numerous and easily observed (see Table 1). Of course the parameters 
are better determined when more than two lines from the same level are used. 



2.4. Fitting the Lines 

In practice however the problem is complicated by the structure of the absorber. 
Most often, an absorption feature is a blend of several components. Indeed 
when the line of sight passes for example through a galactic halo, it can intercept 
several clouds. The relative projected velocities of the clouds are usually not 
very large (typically AV < 150 km s _1 ) compared to the width of the lines 
(5 to 50 km s _1 or larger) and the number of clouds can be large (more than 
ten) . Column densities and Doppler parameters for each component are derived 
from line profile fitting, using the maximum of transitions. This is illustrated on 
Fig. [| for H 1. The fit of a single absorption feature with several components 
is most often not unique; the different solutions can yield column densities 
differing by an order of magnitude. When several transitions from the same 



QSO Absorption Line Systems 9 



Table 1. A few strong atomic transitions 



Ion 


Ao 

(A) 




f 


log(A f) 


log(ASf) 


VI 


1031.927 


0. 


130 


2.128 


5.141 


vi 


1037.616 


0. 


.0648 


1.828 


4.844 


H i 


1215.670 


0. 


.4162 


2.704 


5.789 


i 


1302.169 


0. 


.0486 


1.801 


4.916 


C ii 


1334.532 


0. 


.118 


2.197 


5.323 


Si iv 


1393.755 


0. 


,528 


2.867 


6.011 


Si iv 


1402.770 


0. 


.262 


2.565 


5.712 


C iv 


1548.202 


0. 


.194 


2.448 


5.667 


C iv 


1550.774 


0. 


.097 


2.177 


5.368 


Mg ii 


2796.352 


0. 


.592 


3.219 


6.666 


Mg ii 


2803.531 


0. 


.295 


2.918 


6.365 



ion are used, as the transitions in the Lyman series for H I, better significance 
is achieved. Stronger constraints can be put on models by using several ions. 
Indeed metal lines are narrower than H I lines and column densities are usually 
smaller for metals than for H I, so that metal absorption features are relatively 
less sensitive to blending effects (see Fig. |J). However it must be assumed that 
the decomposition in components is the same for the different ions. Hence, it 
is usually a good assumption to consider that H I, O I, Si II, C II, Al II and 
Fe II on the one hand, and C IV and Si IV on the other arise in the same phase 
and can be fitted together (see Fig. ||). 

Several procedures using \ 2 minimization techniques are available to perform 
Voigt profile fitting of the absorption features. The most popular are FITLY- 
MAN implemented in MIDAS the ESO image processing package (Fontana & 
Ballester f27j) and VPFIT (Carswell et al. §). To visualize the lines XVOIGT 
is also very convenient (Mar |5l]]). Most of the models overplotted on the figures 
have been created with the latter procedure. 



3. THE LYa FOREST 



The numerous H I lines seen in the Lya forest correspond to the imprint left 
by the intergalactic medium in the QSO spectrum. The whole redshift interval 
from zero to the highest QSO emission redshift has been observed. A review 
of the recent works and ideas on this subject can be found in [Q. In the 
following, a few important observations and interpretations are highlighted. 
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Fig. 4. — Upper panel: The fit of a single absorption line (here H IA1215; from 
Fernandez-Soto et al. |^|) is by no means unique; to constrain the model, several 
transitions are used, either from the same ion (in the lowerpanel all the H I transitions 
in the Lyman series are used; from Buries & Tytler [M; note that the model is 
overplotted to the data) or from different ions assumed to arise in the same phase. 
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Fig. 5. — Absorption lines in the z a bs = 3.1825 system toward PKS 2000-330 plotted 
on a relative velocity scale. Parts of the spectrum have been removed when absorption 
lines from other systems were present. It is apparent that (i) the metal lines are much 
narrower than the H I line; note however that metallicity in the system is quite low, of 
the order of 1/100 of solar; (ii) the number of components needed to fit the H I line is 
much larger than for the metal lines; (iii) the position of three metal components have 
been marked by dashed vertical lines; it can be seen that the velocity coincidence is 
very good between Si II, Al n and Fe II. 



3.1. Evolution of the Number Density of Lines 

The redshift evolution of the number of Lya lines per unit redshift is usually 
fitted by a power law, 

N(z) = N (l + z)i (12) 

Note that N depends on the equivalent width; weak lines arc more numerous 
than strong lines (see next Section) . The observed exponent 7 can be compared 
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Fig. 6. — Logarithm of the number density of Lya lines versus redshift from Jannuzi 
et al. |B8|. The break at z ~ 1.5 is apparent. 



with that expected if the forest is due to a non-evolving population: 

7 = 2-i[30 o (l + ^) 3 -2^ o (l + z) 2 ]x[r! o (l + z) 3 -iv o (l + z) 2 + -4^]- 1 (13) 
2 3H 5 

with standard notations and with K Q = o +3^t _ 1 (see [^8]). For a model 
without a cosmological constant, 

7= (l + q z-q )(l + 2q z)- 1 (14) 

In the absence of evolution, the exponent is equal to 1.0 for q D = and 0.5 for 
q = 0.5 at z = 2. Observations show that the Lya forest evolves strongly at 
high redshift (|Q; 7 ~ 2.5 at z ~ 2); that weak lines evolve more slowly than 
stronger one; and that there is a break at z ~ 1.5 in the cosmological evolution 
of the number density of lines (see Fig. ^, pq] ). 

3.2. Column Density Distribution 

The column density distribution for all absorption systems (the Lya forest 
corresponds only to the smallest column densities) at z ~ 2.8 is shown on 
Fig. The best fit for the forest is a power- law d 2 n/diVd z oc N 1,5 . There is 




a deficit of lines at N ~ 10 15 cm -2 and a flatening of the distribution for the 
damped systems [3(| . These features have been interpreted as a consequence 
of the structure of the intergalactic medium and of different ionization effects 
H, [||. It has been shown from this distribution that the Lya forest at 
high redshift contains most of the baryons in the Universe ||66f , a result that has 
been confirmed by other methods |7(], p^j. This suggests that at high redshift, 
the intergalactic medium is the reservoir of baryons for galaxy formation. 



3.3. Physical Conditions 

The density of the gas must be very small. Indeed from coincident absorptions 
along adjacent lines of sight (see below, e.g. [2l[]), it can be shown that the 
gas in the Lya forest is fairly uniform on scales of ~1 kpc. This implies that 
for a typical column density of AT(H i) = 10 14 cm~ 2 , the H I density is smaller 
than 3xl0 -8 cm~ 3 . Even for an ionization correction factor corresponding 
to highly ionized gas, n(HII)/n(HI) = 10 4 , the total density is smaller than 
3xlCT 4 cm- 3 . 

The intergalactic gas is exposed to the effect of the UV-background which is 
the radiation field produced by QSOs and early galaxies at high redshift. It 
is possible to estimate the flux at the Lyman limit using the observed number 
of sources (e.g. [72)), J912 > 2.3xl0~ 22 erg/s/cm 2 /Hz/sr at z — 2.5 if 



the shape of the ionizing spectrum is taken as in 31 . With the above orders 
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of magnitude the ionization parameter^] is U = 0.02 and the gas is indeed 
completely ionized. 

The presence of metals in the Lya forest is an important clue towards under- 
standing how and where the first objects in the universe formed. It has been 
known for long that the amount of metals is small. However given the small 
column densities involved, the question could not be satisfactorily answered be- 
fore the advent of the lOm-class telescopes. Using Keck data, it has been shown 
recently that C IV absorption is seen in all the systems with log N(H i) > 15 
and in half of the systems with log N(H i) > 14.3 @ ||, |l| . A lower limit 
on the metallicity, Z > 10 -3 Zq (H q = 65, q = 0.02), is derived after difficult 
ionization corrections pC| ]. The production of this amount of metals would 
induce enough photons to completely pre-ionized the IGM [Q. However, it 
is much lower than the predictions of metal production by a first generation 
of Type II supernovae p5| . The presence of metals in systems with smaller 
column densities is still controversial (jl7| |50| . 

Another way to estimate the ionizing UV-background is to use the proximity 
effect which is the observed decrease of the number of absorption lines at a 
given equivalent width limit in the vicinity of the QSO. This is interpreted 
as a consequence of the enhanced ionizing flux due to the emission by the 
quasar. The radius of influence of the quasar emission can thus be derived 
from observations. This is the distance from the QSO at which the depleted 
number density of lines reaches to the background value. At this point, to a first 
approximation, the flux from the background is comparable to the flux from the 
quasar Q . However this determination assumes that (i) the number density of 
clouds can be extrapolated toward the quasar (this may be wrong if the density 
of clouds increases in the potential wells of the QSO), (ii) the luminosity of the 
QSO is known from direct observation (questionable if variability occurs), (iii) 
the redshift of the quasar is known (discrepancies between different lines are 
the rule, [p3[). From this it is easy to understand that the typical value of 
J912 ~ 10~ 21 erg/s/cm 2 /Hz/sr derived from these studies is highly uncertain. 
The uncertainties are discussed in ^9) and in more detail in jl5| . 



3.4. Clustering Properties 

The two-point correlation function £ p] has been studied by several authors 
with contradictory conclusions, probably because the number of lines of sight 



was not large enough. Cristiani et al. 18 use 1600 lines and show that there 



is a strong signal on scales Av < 200 km s 1 which increases for increasing col- 
umn densities. However, it is difficult to disentangle the large scale clustering 



(*) Ratio of the density of ionizing photons to the total hydrogen density; since 
riHi x aF/hu a tie x a, with a the ionization cross-section, a the recombination 
coefficient, nm.Hii.e the HI, HII and electronic densities and F the ionizing flux, then 
nmi/nm <x a/a x clJ 

( 2 ) The probability to find a Lya cloud in a volume dV at a distance r from another 
cloud is given by dP — 4>(z)dV[l + £(r)], where $ is the spatial density of clouds. 
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signal from the effects of the internal structure in an absorber |26| . Indeed, the 
column density through a system increases because the number of subcompo- 
nents increases; each component beeing the signature of an individual cloud in 
an extended halo |H| . This is why observers have tried to measure directly the 
dimensions of absorbing clouds. For this, they use pairs or groups of quasars 
with small angular separations on the sky. Correlation of absorption features 
along the different lines of sight indicates that the characteristic size of the 
absorbers is larger than the projected separation of the lines of sight. In the 
spectra of multiple images of lensed quasars with separations of the order of a 
few arcsec J^, [$7|, the Lya forests appear nearly identical, implying that the 
absorbing objects have large sizes (>50 kpc). For larger separations between 
lines of sight, the correlation between the two forests decreases; however after 
random coincidences are removed, an excess of common redshift absorption ex- 
ists even for separations as large as 500 h^ 1 kpc |^0). This suggests very large 
dimensions or at least that the clouds are closely correlated on these scales. 
The data available up to now have been critically analyzed by D'Odorico et 
al. IH who concluded that the absorbers have typical sizes of the order of 
R ~ 350±100/i _1 kpc. The number of suitable groups of quasars is too small 
at present to constrain the structure of the Lya complexes. However, it is 
easy to see that the method is very powerful to study the correlation of bary- 
onic matter at high redshift. The number of observable pairs should increase 
dramatically with the advent of lOm-class telescopes. It is also conceivable, 
although ambitious, to probe the 3D-distribution of the gas in a small field by 
increasing the density of lines of sight. To this aim, quasars with magnitude 22 
should be observed at intermediate resolution. The large amount of telescope 
time needed requires multi-object spectroscopy capabilities on lOm-class tele- 
scopes If a hundred lines of sight could be observed in one square degree, 
the probability distribution function of the matter density could be investigated 
at high redshift. This would be a unique way to study how galaxy formation 
is related to the distribution and dynamics of the underlying matter field (see 
[5S| , [l9| and next Section). The method is illustrated in Fig. |j. 



3.5. Simulations 

Most models of the Lya forest assume that the intergalactic medium is popu- 
lated by diffuse clouds, either gravitationally or pressure confined, and possibly 
embedded in a pervasive medium (e.g. || [76|, |74|, |3(|). A more attractive pic- 
ture arised recently from large-scale dark-matter A^-body simulations including 
a description of the baryonic component: either full hydrodynamic calculations 



]1C| , |34| , [54| , 101 1 limited however to z ~ 2, or an approximate treatment that 
can follow evolution down to z = [5^, [7^, [|. The simulations show that, to 
a very good approximation, the properties of the Lya forest can be understood 
if the gas traces the gravitational potential of the dark matter. In this picture, 
part of the gas is located inside filaments where star formation can occur very 
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z=0.9512 




Fig. 8. — Observation of a quasar pair [^HJ: Q1026-0025A (z cm = 1.438; situated in 
the bottom right corner of the image) and B (z em = 1.526; 36 arcsec or 320^5 1 kpc 
away from A in the north-east direction). North is at the top, east to the left. Lya 
absorptions are seen at z a b s = 0.9486 and 0.9536 in the HST spectrum of A and at 
Zabs = 0.9486 in B. A galaxy is detected with the NTT (ESO-La Silla) at redshift 
z a m = 0.9512 at a projected distance of ~ 450 kpc from the two absorbers. Using 
lOm-class telescopes, it will be possible to observe m ~ 22 magnitude quasar at 
intermediate resolution. The aim is to increase the density of lines of sight in order 
to study the 3D spatial distribution of the gas and to correlate this distribution with 
the galaxies in the field. 
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early in small haloes that subsequently merge to build up galaxies. It is thus not 
surprizing to observe metal lines from this gas. The remaining part of the gas 
either is loosely associated with the filaments and has N(H i) > 10 14 cm -2 , or 
is located in the underdense regions and has N(H i) < 10 14 cm -2 . Much of the 
Lya forest arises in this gas and it is still to be demonstrated that it contains 
metals |l7], ^0|. To clarify these issues, it is important to follow the evolution 
of the Lya gas over a wide redshift interval down to the current epoch (z = 0), 
distinguishing between gas closely associated with the filamentary structures 
of the dark matter and the dilute gas mainly located in the underdense regions. 
The simulations are successful at reproducing the global properties of the Lya 
forest and its evolution (see Fig. ||). They have favored the emergence of a 
consistent picture in which not only the intergalactic medium is the baryonic 
reservoir for galaxy formation, but also galaxy formation strongly influences 
the evolution of the IGM through metal enrichment and ionizing radiation 
emission. More information on the 3D spatial distribution of the gas should be 
gathered from observation of multiple lines of sight to add constraints to these 
models (see previous Section). 



4. METAL LINE SYSTEMS AND GALAXIES 

Detailed studies of the kinematics and metallicities of metal line systems (see 
e.g. [g4| |(3l], ^]) are unique for understanding the structure and evolution of 
galactic haloes. Indeed the connection between apparently normal galaxies and 
strong metal absorption line systems has been explored in some detail recently 
|@@. The fields of quasars known to have metal absorption line systems 
at z a bs are searched for galaxies at redshift z cm ~ z a b s - Up to z ~ 1, there is 
always such a galaxy within ~ 10 arcsec from the line of sight to the quasar. 
At intermediate redshift (z ~ 0.7) Mg n systems with equivalent width w T > 
0.3 A are primarily associated with large R ~ 35ft- -1 kpc haloes of a pop- 
ulation of luminous (~Lg) galaxies which appear not to evolve significantly 
over the observed redshift range. The origin, structure and evolution of these 
large galactic haloes play an important role for the formation and evolution 
of galaxies, as there is a close interaction between the disk where most of the 
stars form and the halo that acts as a reservoir of gas. 

The most striking result of these studies is that the only interlopers (galaxies 
that are close to the line of sight but do not induce absorption in the QSO 
spectrum) have Lk < 0.05 L^. At intermediate redshift, the most important 
criterion for selecting absorbing galaxies seems therefore to be mass rather 
than star-formation rate. However, it has been shown that the spectral energy 
distribution of absorbing galaxies is bluer at higher redshift and that the mean 
equivalent width of [OII]A3727 increases by 40% between z ~ 0.4 and z ~ 0.9 
|3Cfl . This suggests that star- formation activity is also an important factor for 
the presence of large haloes. In any case it seems that the structure of haloes 
is highly perturbed. In particular, it has been shown that the kinematics of 
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Fig. 9. — Left handside panels: 25x25x2 Mpc 3 slice through the simulation box at 
z = 3.28, 1.6 and 0.8; contours delineate regions where JV(H i) > 10 14 cm" 2 . Right 
handside panels: Evolution of the ionizing UV-background flux (calculated from the 
amount of gas cooling below 5000 K in the simulation); H I column density distribution 
at 2 ~ 2; number density of lines with N(H i) > 10 12 , 10 13 and 10 14 cm" 2 versus 
redshift (observational points are overplotted as squares; note the low redshift points 
from ted]). The simulation is from Riediger et al. |75|. 
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Fig. 10. — The Mg II absorption line profiles for ten systems are plotted versus the 
impact parameter of the field-galaxy closest to the line of sight at the same redshift 
as the absorber. It is apparent that there is no correlation between the absorption 
profile and the impact parameter. Although this comparison is very crude since the 
absorbing galaxies have neither the same luminosity nor the same star-formation rate, 
this nonetheless sug gests that galactic haloes have highly perturbed structures (from 
Churchill et al. [p_3[). 



the absorption systems is not correlated with the impact parameter (|13j and 
Fig. |l^). Simple disk or halo models are rejected Q |. 

It is also still unclear which galaxies are associated with weak (w r (Mg n) < 0.3 A) 
metal absorption systems. The characteristic radius of haloes can be estimated 
from the incidence rate of absorption, assuming that every galaxy is surrounded 
by a spherical halo with filling factor unity, and that the relation between lumi- 
nosity and halo radius derived from Mg n galaxy detection can be extrapolated 
to small equivalent width. The detection limit for weak metal lines has been 
pushed to very low values with lOm-class telescopes. Mg n systems are ex- 
pected to trace the gas with highest density. The number of Mg II systems 
per unit redshift at 0.4 < z < 1.4 with W t > 0.02 A is dN/dz = 1.74±0.10 
|14}| . This is four times larger than the number of Lyman limit systems and 
corresponds to 5% of the number of Lya forest lines with W T > 0.1 A. This 
implies a halo radius of the order of 63ft- -1 kpc for L* galaxies. This suggests 
that either normal galaxies are surrounded by huge haloes or that weak systems 
are associated with a population of dwarf galaxies. 
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5. DAMPED SYSTEMS 



Damped Lya (hereafter DLA) systems are defined as systems with hydrogen 
column density N(R i) > 2xl0 20 cm~ 2 . This definition is artificial since 
damped wings appear for lower column densities (iV(H i) > 10 19 cm -2 ; see 
Section 2.3 and Fig. ||). It has been introduced assuming that these lines 
should be characteristic of galactic disks at high redshift |)7J . Another reason 
is that damped Lya surveys were performed at low resolution. The equivalent 
width w ohs (z ~ 2.5) > 17.5 A for iV(H i) > 10 20 cm" 2 . The probability that 
such a strong absorption feature is the result of blending is small (see Table 3 
of ]97]]). It is clear that this definition may introduce a systematic bias in the 
discussion of what is the nature of these systems and it would be most valu- 
able to compare the properties of systems with 19 < log A(H i) < 20.3 and 
log N(H i) > 20.3. 

For log JV(H i) > 19, the optical depth at the Lyman limit is large enough so 
that hydrogen is neutral. The gas is either cold (T < 1000 K) and contains 
molecules (e.g. |83|| ) or warm (T ~ 10 4 K) for the highest or lowest column 
densities respectively |33| . As a consequence of the shape of the column density 
distribution, d 2 n/dNdz oc A - * 3 with (3 ~ 1.5, most of the mass is in the systems 
of highest column densities. The number density of the latter decreases with 
time presumably as a consequence of star-formation ( |97|, [l2) ; see however (9l) ) . 
Indeed, the cosmic density of neutral hydrogen in DLA absorbers at z ~ 3 is 
similar to that of stars at the present time |89| . 

Metallicities and dust content have been derived from zinc and chromium ob- 
servations |(37], ^38). The [Zn/Cr] ratio is an indicator of the presence of dust if it 
is assumed, that, as in our Galaxy, zinc traces the gaseous abundances whereas 
chromium is heavily depleted into dust grains (see jijj and |59| for a critical 
discussion of this assumption). The typical dust-to-gas ratio determined this 
way is about 1/30 of that in the Milky- Way |68[ |95|. This amount of dust could 
bias the observed number density of DLA systems fjj . Metallicities are of 
the order of a tenth solar and tend to decrease from z~2toz>3(|69 ) 7j and 
Fig. |lTJ) . At any redshift however, the scatter is large and it may be hazardous 
to draw premature conclusions from the small samples available. 
Recently, Prochaska & Wolfe [[n| have used Keck spectra of 17 DLA absorbers 
to investigate the kinematics of the neutral gas from unsaturated low-excitation 
transitions such as Si 11AI8O8. They show that the absorption profiles are in- 
consistent with models of galactic haloes with random motions, spherically 
infalling gas and slowly rotating hot disks. The CDM model Q is rejected 
as it produces disks with rotation velocities too small to account for the large 
observed velocity broadening of the absorption lines. Models of thick disks 
(h ~ 0.3i?, where h is the vertical scale and R the radius) with large rotational 
velocity (v ~ 225 km s _1 ) can reproduce the data. In a subsequent paper how- 
ever, Haehnelt et al. f32|| use hydrodynamic simulations in the framework of a 
standard CDM cosmogony to demonstrate that the absorption profiles can be 
reproduced by a mixture of rotational and random motions in merging proto- 
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Fig. 11. — Zinc metallicity (relative to solar) in damped systems versus redshift (from 
Pettini et al. (6§]). 

galactic clumps. The typical virial velocity of the haloes is about 100 km s _1 . 

The issue of whether DLA systems originate in thick disks has been ques- 
tioned previously. In particular, the metallicity distribution of DLA systems is 
incon sistent with that of stars in the thick disk of our Galaxy ( |p9|| ; see however 
|l00| ). Arguments in favor of DLA systems being associated with dwarf galax- 
ies have been reviewed by Vladilo |95f| . However, it has been shown recently 
that DLA systems at intermediate redshift are associated with galaxies of very 
different morphologies Q] . This strongly suggests that the objects associated 
with high-redshift DLA absorbers are progenitors of present-day galaxies of all 
kinds (see @). 

6. CONCLUSION 

The amount of information derived from studies of QSO absorption line systems 
has increased tremendously during the last few years. This largely results from 
the improvements in instrumental sensitivity, and spectral resolution and the 
extension of the accessible wavelength range. Coupled with the development of 
sophisticated A^-body simulations, this progress has favored the emergence of 
a consistent picture in which not only the intergalactic medium is the baryonic 
reservoir for galaxy formation but also galaxy formation strongly influences 
the evolution of the IGM through metal enrichment and ionizing radiation 
emission. 

One of the most exciting projects for the next few years may be to probe the 3D 
spatial distribution of the gas at high redshift. A number of large surveys are 
underway or planed, which will provide a large number of bright (my < 19) 
quasars for investigation of large scales. A few small fields will be searched 
for deeper samples of quasars (my < 23). If a hundred lines of sight can 
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be observed in one square degree, the probability distribution function of the 
matter density could be investigated at high redshift. This would be a unique 
way to study how galaxy formation is related to the distribution and dynamics 
of the underlying matter field. 
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